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Abstract 

This paper presents a numerical study about the performance of a Beta Stirling solar thermal engine system. This 
system is composed of a solar collector box connected to a regenerator hydraulic system and a transmitting power sys¬ 
tem. The objective of the system is to offer a new alternative to help solving stagnant water pollution in hot countries 
like Thailand by circulating water in canals, lakes, ponds etc. for aeration using solar energy. 

The purpose of this study is to determine the power output and actual heat transfer on the performance of the solar 
thermal engine. The solar thermal engine is analyzed using a mathematical model based on the first law of thermody¬ 
namics for processes with finite speed, with particular attention to the energy balance at the receiver. The result of cal¬ 
culations showed that the regenerator volume and phase angle must be chosen carefully to fulfill the requirement that 
total fluid mass in the system is constant and to obtain maximum power output throughout the day. 

© 2004 Elsevier Ltd. All rights reserved. 
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1. Introduction 

Environmental impact and depletion of mineral re¬ 
sources such as gas oil and lignite are prompting a re- 
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examination of an alternative to these resources. 
Renewable energy sources seem to be a good alternative. 
The power needed to drive turbine to circulate water of 
canals, lakes ponds, etc. must come from local energy 
sources. Besides electrical motor power, solar and wind 
energy, and to some extent biogas, may provide the solu¬ 
tion to this problem. 

Solar energy is the most fascinating and promising 
renewable energy source. Thailand is a country situated 
within an equatorial belt having hot and humid climate. 
The average daily intensity of solar radiation is 17.5MJ/ 
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Nomenclature 



Ac 

collector area (m ) 

Opiate 

mean plate temperature (K) 

dx 

phase angle 

71 

intermediate low temperature sink (K) 

n 

engine speeds (rps) 

U\_ 

coefficient of the thermal energy lost (W/ 

p 

A mean 

mean engine pressure (Pa) 


m 2 K) 

P 

1 min 

minimum engine pressure (Pa) 

V 

volume (m 3 ) 

P 

* max 

maximum engine pressure (Pa) 

fsE 

swept volume of expansion piston or dis- 

PWi 

theoretical output engine power (W) 


placer piston (m 3 ) 

Qh 

heat transfer (hot side) (kJ) 

He 

swept volume of compression piston or 

Qc 

heat transfer (cold side) (kJ) 


power piston (m ) 

R 

gas constant 

V B 

overlap volume (m ) 

S 

entropy of the working fluid (kJ/kgK) 

Udc 

compression dead volume (m ) 

s 

solar radiation absorbed by a collector (W/ 

VjDE 

dead volume of expansion space (m ) 


m 2 ) 

Vr 

regenerator volume (m ) 

t 

temperature ratio of T E /T C (K/K) 

X 

crankshaft angle (degree) 

r, 

ambient temperature (K) 

W c 

energy of the compassion (J) 

Tc 

constant cold temperature of the working 

W E 

energy of the expansion (J) 


fluid in the compression space (K) 

W\ 

energy output (J) 

Te 

constant hot temperature of the working 
fluid in the expansion space (K) 

Oe 

thermal efficiency 


m day (Namprakai et al., 1989). Therefore, it is neces¬ 
sary to consider the use of this freely available energy 
in various fields such as the conservation of food prod¬ 
ucts, the habitat, and its transformation into mechanical 
energy. In recent decades, some solar pump and solar 
thermal engines operating on the principle of thermody¬ 
namic conversion scheme have been built and tested 
extensively throughout the world (Kaushik and Kumar, 
2000). These systems used the thermal energy produced 
by the solar energy to drive a conventional engine for 
pumping water (Tallbert et al., 1978). However, many 
new designs have been introduced recently and the 
developments in this field need to be updated. The aim 
of this paper is to study the performance of a solar ther¬ 
mal system for circulating water in canals, lakes ponds 
etc. for aeration. In fact, water pollution, is not only 
due to garbage and waste disposed into the water, but 
also as a result of poor circulation of water or water 
stagnation in many canals, lakes and ponds. This water, 
usually in urban and suburban areas, is enclosed and has 
weak wind speed to help in the circulation. 

2. Literature review 

The objective of this review is to provide a basic 
background and some review of existing literature on 
solar-powered Stirling engine technology that uses low 
temperature technology. A number of Stirling engine 
configurations and designs, including the engine’s devel¬ 
opment, are provided and discussed. 


The extent of research on solar Stirling-cycle for 
power production has been limited to fraction of horse¬ 
power (Farber and Prescott, 1965). The hot-air engine 
developed by Farber and Prescott in 1965 focused on 
area where air was heated and its expansion pushes the 
piston down. In the down-stroke of the piston, the dis¬ 
placer moved to the left by a linkage. On the up-stroke 
of the piston, the displacer moved to the right and all 
the hot air was at the left section of the cylinder and lost 
heat to the cooling water. An efficiency of about 9% was 
obtained at lOOrpm with a break horsepower of about 
0.2. And after that Beale et al. (Beale et al., 1971) in 
1971, considered a free cylinder containing a heavy pis¬ 
ton which remained essentially stationary and a light dis¬ 
placer which moved under the influence of pressure 
differential between the work space and the bound space 
in the cylinder. The cylinder moved under the influence 
of the same pressure differential and performs work 
against an external load. Since the entire pressure enclo¬ 
sure moved as a unit, the free cylinder engine can deliver 
work from a completely sealed working gas and since 
there are no internal bearing loads, or mechanical link¬ 
age or gears, etc., no lubrication was needed. The engine 
should be inexpensive and should have a long life. A 
Fresnel lens focused the solar radiation on the hot space 
of the engine and a double acting water pump utilized the 
developed mechanical energy. An efficiency of about 9% 
was obtained at lOOrpm with a power of about 149 W. 

Stirling engine theoretically could achieve the highest 
possible energy conversion efficiency of all heat/power- 
engines. Using practically any kind of energy, they can 
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be employed with extremely low emissions as generator, 
motor, heat pump or cooling system. In spite of all these 
advantages, Stirling engines still have not made the mar¬ 
ket breakthrough as a mass product. A general analysis 
of finite time thermodynamics of a Stirling heat engine 
was presented with finite heat capacity of external reser¬ 
voirs, regenerative losses and finite effectiveness of each 
of the heat exchangers (i.e. a high temperature heat ex¬ 
changer and a low temperature heat exchanger) (Fette, 
1995). They obtained the expressions for maximum 
power output and computed the corresponding thermal 
efficiency. The effect of operating temperatures, the effec¬ 
tiveness of the regenerative heat exchanger on the heat 
transfer (g H and g L ) to and from the heat engine, the 
regenerative heat transfer (g R ), the maximum power 
output ( P ) and the corresponding thermal efficiency of 
the cycle have all been studied. They considered the inlet 
source temperature and sink temperature as 1300K and 
300 K. They also used effectiveness of each heat exchan¬ 
ger and the regenerator in the range from 0.4 to 1, and 
the capacitance rates of the heat source/sink reservoirs 
in the range from 0.30 to 1.80kW/K. Their analysis 
showed that power output and thermal efficiency were 
62.55kW and 43.94%, respectively. 

3. Description of solar thermal engine 

The solar thermal engine is essentially composed of a 
solar collector, a displacer, a cooling system, a regenera¬ 
tor, a hydraulic and a transmitting power system (Fig. 1). 

The solar collector is a box, where the absorber is an 
aluminum plate, and the cover is a transparent acrylic 


(2-layer fluor-acrylic) polymer. The thermal radiation 
absorbed by the absorber is used to heat the liquid ethyl 
fluorocarbons (R-ll) and the absorbed energy set the 
fluid in motion by thermo-siphon action and the liquid 
is transformed into vapor. We are aware that R-ll is 
not environment friendly. However, due to its proper¬ 
ties, which are suitable to operate the engine, we conduct 
our theoretical study using R-ll so that we could formu¬ 
late a basis for the design of such a thermal engine. The 
outside of the box next to the regenerator is covered with 
insulation material to prevent heat loss to the outside. 

The displacer moves up and down part of the box. 
When the displacer moves up it pushes the working 
vapor against the inside of the hot absorber plate. A dis¬ 
placer membrane, made of a flexible and heat-insulation 
fabric and attached to the outside ring of the displacer, 
connects the displacer to the interior lining on the out¬ 
side of the solar collector box. 

The cooling system located below the regenerator 
further cools down the temperature of the working fluid 
before it enters the lower and cool part of the solar col¬ 
lector box. The cooling of water-fed pipes absorbs heat 
from the working fluid flowing through. Cold water is 
circulated though the cooling system and thus keeps this 
area cool. 

The hydraulic system is a heat exchanger between 
water and ethyl fluorocarbon vapor. The transmitting 
power system transforms the translation movement into 
a rotating movement. The function of the power piston 
is to convert the expansion of the working fluid at high 
pressure and to compress the working fluid at low tem¬ 
perature and to transfer this conversion into motion by 
means of a crankshaft and flywheel. 


Qext 


Transparent 



Fig. 1. Schematic diagram of the solar thermal engine. 








P. Kerdchang et al. / Solar Energy 78 (2005) 518-527 


521 


When the displacer has reached the bottom position, 
cool vapor of ethyl fluorocarbons fills the box. So, the 
temperature and the pressure of the vapor reached its 
low value (about 20-30 °C). When the displacer moves 
up, the process is reversed. Cool vapor of ethyl fluoro¬ 
carbons is pushed through the regenerator picking up 
the heat stored there flowing along the inner hot side 
of the box, and its temperature increases. As soon as 
the displacer has reached its highest temperature (about 
70-100 °C) there is a fluctuation of the pressure in the 
box. The bottom of the box is attached to the absorber 
structure through a flexible membrane coupling system. 
Therefore, it can serve as a large-size power piston in re¬ 
sponse to the cyclical pressure fluctuation inside the box. 
The pressure fluctuation is used to generate mechanical 
work in the following way. The displacer moves the 
vapor cyclically from the hot to the cool section of the 
box. So, the evolution of the pressure flow is cyclical 
and subsequently the displacer movement produces a 
mechanical work. For half a cycle, the regenerator stores 
heat of the working fluid vapor temporarily and the 
cooling system is used to cool the vapor before it flows 
in the lower part of the box. 



Fig. 3. T-S diagram of the solar thermal engine cycle. 


4. Theoretical model 

4.1. The thermal engine 

A simplified model (Hsu et al., 2003) as shown in Fig. 
2 can represent the heat process of the solar thermal en¬ 
gine. The T-S diagram (Blank and Wu, 1995; Berrin 
Erbay and Yavuz, 1997) given in Fig. 3, shows the effects 
of heat transfer, regeneration time, and imperfect regen¬ 
eration on the performance of the irreversible Stirling 
engine cycle. 


The hypotheses for modeling the thermal engine are 

common. They are as follows: 

• The expansion and the compression processes are 
isothermal. 

• The vapor is considered as an ideal gas. 

• There is a perfect regeneration. 

• The closed system has constant fluid mass. 

• The regenerator gas temperature is an average of the 
expansion gas temperature (T E ) and of the compres¬ 
sion gas temperature ( T c ). 



Fig. 2. The simplified heat transfer model of the solar thermal engine. 
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• The expansion space (F E ) and the compression space 
(V c ) change according to a sine curve. 


The variation of the volume ratio in the engine sys¬ 
tem is performed by varying the volume of working liq¬ 
uid within the collector box below the working fluid 
vapor cushion. 

The thermal efficiency of the process is defined as the 
ratio between 44 IF” the external work output done per 
cycle and “0 ext ” the external heat supplied by solar en¬ 
ergy per cycle to the engine (Costea and Feidt, 1998). 


r] 


W 

Qqx t 



Concerning all phases the efficiency of the process of 
Fig. 2 is 


0 


Q\ 4 fil-2 

03 4 — (62-3 + 24-l) 



In this ideal Stirling process, the regenerator effi¬ 
ciency is assumed to be equal to 100%, where 
0 ext = external heat supplied by cycle (kJ); Q x _ 2 = heat 
rejected by the cooler at the compression phase at 
44 T c ” (kJ); 02-3 = isochoric heat supplied by V mm (kJ); 
0 3 ^ = external heat supplied by the expansion phase 
at 44 T E ” (kJ); 04^i = isochoric heat rejected by F max (kJ). 

So 02_3 + 04_i = 0 and the external supplied heat for 
the ideal process is 0 ext = 0 3 ^. So that the ideal effi¬ 
ciency is 




63-4 Q 1-2 

63-4 



Following the analysis conducted by Hirata (Hirata, 
1997) on Beta Stirling engine, the momental expansion 
volume (F e ) and compression volume (F c ), described 
at a given crank angle, based on above assumptions are: 

Fe = —(1 — cos(x)) + Fde (4) 

Vc = fp(l -cos(x)) +fp{l -cos(x-dx)} 

+ V DC - K b (5) 

where x = crankshaft angle (degree); V = volume (m ); 
F se = swept volume of expansion piston or displacer 
piston (m ); F S c = swept volume of compression piston 
or power piston (m ); F B = overlap volume (m ); F DC = 
compression dead volume (m ); F DE = dead volume of 

2 2 

expansion space (m ); F R = regenerator volume (m ). 
The overlap volume (F B ) is given by 


The volumes of the expansion and compression cylin¬ 
der at a given crank angle are determined at first. The 
momental volumes are described with a crank angle 
(x). This crank angle is defined as x = 0 when the expan¬ 
sion piston is located at the top most position (top dead 
point). 

The total momental volume (V) is the sum 

V = F e + F c + F r (7) 


The working fluid pressure ( P ) is deduced from the 
expression: 


P 


p Vi — c : 

1 mean v 1 


1 — C • cos(x — a) 
F max ( 1 - C) 

1 — C • cos(x — a) 


F mm ( 1 + C) 

1 — C • cos(x — a) 



The energy of the expansion (fF E ) and compression 
(fF c ) spaces are expressed respectively by 


Fmean FSE^F • sin($) 

1 + Vi - c 2 

P mill KSE^ * 0 • SI 11 ( (7) 

1 + V1 - c 2 * Vi-c 

Pm axK SE 7i • C ■ sin (a) ^ Vl - C 

1 + Vi - c 2 * Vi + c 



w, 


C 


-7 


PdV c = - 




Fmean F SE ' TC ' C • t • sin({2) 

1 + Vi - c 2 


P min F SE 71- C t- sffi(a) \/\+C 

— -* 


l + Vi - c 


VY^c 


FmaxF se ’ tc ‘ C • t • sin(^) "\/l F 

l + Vi - c 2 * vTTf 


( 10 ) 


c = 


B 

~S 


(10a) 


4 tx 

S = t + 2lTe)e + z -b F + 2 Xdc H - 1 

A I t 


(10b) 


B 


P + 2 (t- i)V cos(dx) 




(10c) 


a — tan 1 


F • sin(dx) 
t + cos(dx) + 1 


(lOd) 


where P mean = mean pressure (Pa); P m \ n = minimum 
pressure (Pa); P max = maximum pressure (Pa); t = tem¬ 
perature ratio of T E /T C . 

For one cycle, the energy is equal to 

Wi = W E + Wc ( 11 ) 



where W E = energy of the expansion (J); W c = energy of 
the compassion (J); W\ - energy output (J). 

The engine power (P W[ ) and the thermal efficiency 
(rjo) are defined respectively: 
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W 

P w . = Win; ti e =~y=l -t 

VV e 

where n - engine speeds (rps). 
4.2. The solar collector box 



where f/ L = coefficient of the thermal energy lost (W/ 
m“); S - solar radiation absorbed by a collector (W/ 
m 2 ); T F m = mean plate temperature (K); T a = ambient 
temperature (K); A c = solar collector box area (nr). 


An energy balance that indicates the distribution of 
incident solar energy into useful energy gain, thermal 
losses, and optical losses describes the theoretical con¬ 
sideration of system of the solar collector box. The solar 
radiation absorbed by solar collector box is equal to the 
difference between the incident solar radiation and the 
optical losses. In steady state, the useful energy output 
of a collector is then the difference between the absorbed 
solar radiation and the thermal loss (Duffie and Beck¬ 
man, 1991): 


Q u = MS - u L (r P , m - r d ) 



5. Results and discussion 

To discuss the performance of this solar thermal en¬ 
gine, we use the data for the solar radiation and ambient 
temperature on the date 2/02/2001 in Bangkok, Thailand. 
The parameters used in the simulation were 

V SE = 0.0269 (m 3 ), V sc = 0.0269 (m 3 ), 

F dc = 0.0237 (m 3 ), V DE = 0.0237 (m 3 ), 

V R = 0.05 to 0.5 (m3), 

T e = 50 + 273 (K) to 200 + 273 (K), 

T c = 30 + 273 (K), P mean = 180 (kPa), 

A c = 0.39 (m2). 



Volume (m 3 ) 


Fig. 4. Simulated P-V diagram for phase angle 90° at various 
temperature. 



Volume (m 3 ) 


Fig. 5. Simulated P-V diagram for phase angle 140° at various 
temperature. 



Crank shaft angle (degree) 


T e =50°C 

-T e =75 °C 

- T e =100°C 

• T e =200 °C 

* Volume (m 3 ) 


Fig. 6. Pressure and volume when phase angle is 90° (Fr = 0.5 m ). 
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Fig. 7. Pressure and volume when phase angle is 140° (F R = 0.5m ). 
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Fig. 8. Graphs of error of total mass at different phase angle and various T E : (a) 200°C; (b) 100°C; (c) 75°C; (d) 50°C. 
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Fig. 4 shows the simulated P-V diagram for 90° 
phase angle at various gas expansion temperature, and 
Fig. 5 shows that for 140° phase angle. Figs. 6 and 7 
show the pressure variation versus the crankshaft angle 
at different gas expansion temperature. Obviously, the 
higher is the temperature, the higher the pressure and 
the larger the difference between the compression and 
expansion portions of the curves of the cycle; i.e., higher 
power work output. 

One of the assumptions of Schmidt’s theoretical 
model (Hirata, 1997) was that total mass of gases in 
the system was constant because it was a closed system. 
A check of this assumption by simulation showed that 
total mass would not be the same when the phase angle 
was varied, as shown in Fig. 8. The graphs show that the 
phase angle at which the minimum error in total mass 
varies with T E (hot-side temperature) and independent 
of F r . The values of the angles are shown in Table 1. 
However, V R affects the slope of hot-side temperature 
Vs phase angle curve. When V R < 0.1 m , the slopes 



(a) Phase angle (degree) 



(C) Phase angle (degree) 


Table 1 


Phase angle at which error in total mass of vapor is minimum 
(1.25%) “ 


Hot-side 
temperature (°C) 

Phase angle at which total 
mass is approximately constant 

50 

155 

75 

144 

100 

138 

200 

122 


are steeper than when V R > 0.5 m . The maximum error 
was estimated to be 1.25%. 

Since solar radiation varies throughout the day, the 
hot-side temperature can be expected to vary between 
50 °C and up to about 200 °C. To accommodate this 
temperature variation, phase angle must be varied be¬ 
tween 122° and 155°. It will be impractical to vary the 
phase angle according to the solar time through out 



(b) Phase angle (degree) 



(d) Phase angle (degree) 


Fig. 9. Simulated power outputs at varying phase angle and different T E : (a) 100°C; (b) 200 °C; (c) 50°C; (d) 75 °C. 
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Time (h) 


Fig. 10. The relationship between solar radiation and simulated power output (V R = 0.5 m 
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Fig. 11. Relation between solar radiation and engine speed (F R = 0.5m ). 


the day. A compromise must be sought which could be 
used throughout the day. F R could be chosen so that 
error in total mass will be at an acceptable minimum. 
Fig. 9 shows simulated power output graphs. The graphs 
show that output power is a function of hot-side temper¬ 
ature and F R . Power output is low at high F R , but it is 
high at low F R . Fortunately, the effect of F R on power 
output is low if F r is greater than 0.5 m , but when 
F r is equal to 1, output is only half of that when F R 
is 0.5. Therefore, F R = 0.5 m 3 is suitable for this simula¬ 
tion. Then, the compromised choice of the phase angle 
would be around 140° with the associated error in mass 
of 1.2% in the cyclical process. 

Fig. 10 shows the relationship between solar radia¬ 
tion and simulated power output and efficiency of solar 
thermal engine. It was the result of calculation by using 
the solar radiation data and estimated temperatures (r E ) 
and ( T c ). Temperature (T E ) was estimated as 75% of 


(Tpiate) and temperature (T c — 30) was taken approxi¬ 
mately as 75% of water temperature. 

Fig. 11 shows the relationship between solar radia¬ 
tion and the rpm available from the solar engine consid¬ 
ered. The maximum rpm achievable was 120. Since solar 
radiation is low during morning hours and evening 
hours, the speed of the solar engine would be reduced 
accordingly. 

6. Conclusion 

Although, many new designs have been introduced 
recently and the developments in this field need to be 
updated, a new solar thermal engine for circulating water 
of canals, lakes, fish and shrimp ponds, and irrigation 
reservoirs was studied numerically. The power output 
versus the energy output characteristic, heat conductance 
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characteristic from solar radiation to conceivable heat on 
the collector box, and imperfect regeneration coefficient 
characteristic were simulated. Results showed that the 
power and consequently the crankshaft velocity of this 
system depend highly on the thermal performance of 
the solar collector. To meet the assumption of constancy 
of mass in the system, regenerator volume must be care¬ 
fully chosen, and the phase angle must be set at a value to 
give maximum power output for the whole day. Based on 
simulation results, an experimental prototype of this en¬ 
gine is under deployment, and results will be published in 
a coming paper. This system once validated with experi¬ 
ments, is expected to help reducing stagnant water pollu¬ 
tion in both rural and urban areas and to promote the use 
of solar energy as well. 
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